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ON THE FORM AND POSITION OF THE SEA-LEVEL AS DEPENDENT ON 
SUPERFICIAL MASSES SYMMETRICALLY DISPOSED WITH RE- 
SPECT TO A RADIUS OF THE EARTH'S SURFACE. 

By Mr. R. S. Woodward, Washington, D. C. 

1. A considerable class of problems concerning the sea-level is that in which 
the attracting or disturbing mass is symmetrically disposed about a radius of the 
earth's surface and is situated on or near that surface. Examples of such masses, 
proximate at least, are portions of the earth's crust, like a continental table-land, 
and the ice-caps which are supposed to have been of great extent and thickness 
during certain epochs of the earth's history. 

It is proposed in the following paper to develop the theory of the solution 
of this class of problems so far as is necessary to render practicable the numeri- 
cal evaluation of the characteristic effects of the disturbing mass in any special 
case. 

2. The solution of all problems of the kind outlined above depends on the 
principle of hydrostatics that the potential of the forces producing a free liquid 
surface in equilibrium has a constant value for all points of that surface. In the 
case of the earth, if the potential of all the attractive forces acting on a unit mass 
at any point of the sea-surface be denoted by P, the distance of the point from 
the axis of rotation by / and the velocity of rotation by w, the form of the sur- 
face will be completely defined by the equation 

P -f ^/ V = a constant. (i) 

The surface defined by this equation is now commonly called the geoid. It, 
as represented by the ocean surface, is a real surface and does not necessarily 
coincide with the earth's ellipsoid, which is an ideal surface of reference. 

3. The exact value of P in (i) is a complicated function of the densities of 
the element-particles of the earth and of the co-ordinates of those particles and 
the attracted point. For the present purposes, however, it will suffice to con- 
sider P due to a centrobaric sphere of equal mass and volume with the earth 
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and concentric with the earth's centre of gravity. Since we shall only deal with 
relative positions of any point on the sea-surface, the potential due to centrifugal 
force, which is represented by the second term in (i), may be neglected. 

If a and b denote the equatorial and polar semi-axes respectively of the 
earth's ellipsoid, and r^ the radius of the sphere just referred to, 

r„=r^. (2) 

The known values of a and h (Clarke, 1866) give r„ = 20902394 English feet, 
or in round numbers 21000000 feet. The surface of the sphere thus defined may 
be regarded as the surface assumed by a thin film of sea-water covering a nu- 
cleus whose mass plus the mass of the film equals the earth's mass. We shall 
call this ideal surface the undisturbed surface. With respect to it the real sur- 
face of the earth lies partly without and partly within ; but, so far as small rela- 
tive changes in sea-level are concerned, it is practically immaterial whether we 
refer to the actual, closely spheroidal surface or to the simpler spherical one. 

4. Let M^ mass of the earth, 

p^ = mean density of the earth. 
Then J/=fr„V„; (3) 

and the equation to the undisturbed surface is 

— = W^Pm = C^, a constant. (4) 

'0 

Suppose now a new mass m, of density p (positive or negative), placed in 
any fixed position relative to the undisturbed surface. The resulting new sea- 
surface will then differ from that defined by (4). To determine this difference, 
let Fbe the potential of the disturbing mass m at any point of the disturbed 
surface, and let v denote the elevation or depression of this point with respect to 
the undisturbed surface. The equation of the disturbed surface will then be 

M 
— I 1- F = C2, a constant. (5) 

r^ + v 

The difference of this and (4), to terms of the first order inclusive in v, is 
M 

whence, putting ''0 ^ ^2 — ^i> 

we have v = {V-Vf^. (6) 

Since (if we omit the unit of attraction, which disappears by division from 
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our expression) Mjr^^ ^ g, the velocity-increment at the surface of the earth due 
to the earth's attraction, (6) may be written 

. = ^1^^. (60 

5. Vf, in the last two equations is the value of Fwhere z' = o, or the value 
of F along the line of intersection of the disturbed and undisturbed surfaces. If 
we put' 

-+%=^^=f (7) 

This equation represents the elevation of the disturbed surface above a 
spherical surface of equal potential, whose value is 

M 



r. — v 



— 62, 



since the difference between this and (5) gives (7). 

The constant J^ may be determined from the obvious condition that the 
disturbed and undisturbed surfaces contain equal volumes. It should be observed 
also in this connection, that V^ is never of a lower order than V, and cannot 
therefore be neglected in comparison with V. This remark is important since 
some writers* have neglected V^ and arrived at the erroneous equation 

_V 
^~ g 
instead of (6') or (7). 

6. It is evident that the equations derived above, (5) to (7), will hold true if 
the disturbing mass m h& a part of the earth's mass so long as the ratio mjM 
may be neglected in comparison with unity. Thus, in may represent the mass 
of a continent, the deficiency in mass of a lake or lake-basin, etc. 

7. The next step requires the determination of the potential Fof the at- 
tracting mass for any point of the disturbed surface, whether without or within 
the circle which we have assumed to define the boundary of the mass. 

In order to derive an expression for V, let the rectangular and polar co-or- 
dinates of any point of the disturbing mass be defined by the usual relations, viz.: — 

x=^r cos d cos X, 
y ^r cos d sin \, 
z = rsmd\ 

♦Notably Archdeacon Pratt. See his Figure of the Earth, 4th edition, articles 200 and 213. 
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in which d and k correspond to polar distance and longitude respectively, the 
position of the origin being arbitrary. With reference to the same origin, let the 
co-ordinates of the attracted point on the sea-surface be 

x' ■=-r' cos d' cos V , 
y' = r' cos 6' sin X' , 
z' = r' sin 6'. 

If D denote the distance between the attracting and attracted points and 

cos (l> = cose cos d' + sin d sin d' cos {X — X'), (8) 

D"* = r^ -\- r'^ — 2rr' cos <p = {r — r'f + ^rr' sin' ^(p. (9) 

The volume-element of the attracting mass is 

dxdy dz^ r^dr sin d dO dX. 

Hence if p denote the density, supposed uniform, of the attracting mass, a 
general expression for the required potential is 

-r Vr sin d dd dX 



^//J- 



D 



(10) 



We must now evaluate this integral. Taking the centre of the sphere of 
reference as origin of co-ordinates, let 

L^'"''" (II) 

and r' = r^-\-v, ^ ' 

in which u and v are small quantities relative to r^yV being the same of course as 
defined by (6). Premising what will be proved hereafter, namely, that quantities 
of the order «/>•„, vlr^, and (« — vf j r^ may be neglected, equation (9) gives 

Z» = 2r„sin^^ (12) 

From the first of ( 1 1) dr=du, 

and r'=r,\ ^'^' 

to terms of the order u / r. 

Without loss of generality we may assume the line from which d and 6' are 
reckoned to pass through the attracted point, and the plane from which X and X' 
are reckoned to pass through the attracted point and the centre of the attracting 
mass. In this case d' = o, and X' = o, and (8) gives 

4> = d. 
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lOI 



By means of this relation and the values in (12) and (13) equation (10) be- 
comes 



V= r^pC ( (du cos \e dO dX. 



(14) 



Assuming for the present that the attracting mass is of uniform thickness, 
h, the limits of u will be o and h. Let the limits of d, which are obviously func- 
tions of X, be di and 0^. The limits of X are evidently equal in magnitude, but of 
opposite signs. Hence we have 



k e. 



2r^p \ \ \ du cos ^d dO dX ■■ 



= 4^0' 



o 6, o 



A 

kpj^ (sin 1^2 — sin -^i) dl. (15) 



To complete the evaluation of (15) it will be convenient to change variables. 
Consider the spherical triangles formed by the attracted and attracting points, the 
centre of the attracting mass, and the points in which the arc d cuts the circle 
bounding the mass. Thus, in figures i and 2 let P be the attracting and A the 
attracted points, C the centre of the attracting mass and BDG the bounding 
circle. Then, 

d = AP 



and 



X = BAC. 




c 

Fig. I. 

Draw CE perpendicular to AB and put 

AC = a, BC=^, 

PE=s, B£=So, 

CE = p, AE=q. 

From either figure d ^q -\- s, 

^2 = ^1 + ^0> 




whence 



sin ^d^ — sin ^^1 = 2 cos ^q sin ^s^. 



(16) 
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The right-angled triangles of either figure give 

cos a 



co^q- 



cos/' 



The first two of (17) give 



COS 5 , , 

cosi„ = -, (17) 

" cos/ ^ " 

sin / = sin a sin X. 



, 1 , cos a 

2 cos' \q= \ 



cos/' 



, , cos B 

2 cos^ ^ .f„ = I — ^ • 



cos/ ' 
whence 2 cos \q sin ^ = ^K'^"^ ^ + '^^^ ""^^^^^ ^ " ^'^^ ^)] . (18) 

From the last of (17) 

,, cos p dp , , 

"■'•==— n — z / — 2 — X- (19) 

l/(cos^/ — cos^ a) ^ ^' 

Now the last of equations (17) and the diagrams show that the limits of/ cor- 
responding to the limits of X are o and a , or o and /9, according as the attracted 
point is within or without the circle bounding the attracting mass. Hence, if we 
denote the potentials in the two cases by V^ and V^ respectively, the equivalents 
in (15), (16), (18), and (19) give 

o 

«</5 

o 

8. The integrals in these equations are in general elliptics of the third 
species. They may be evaluated by the usual processes applicable to elliptics, by 
mechanical quadrature, or by series. The integral in (20) presents some apparent 
difficulty, since the element-function is infinite at the upper limit except in the 
case a := j8. Again in case a := o, the mtegral assumes this anomalous form 



U \}^^] ■* 
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the value of which is w sin ^/3, as may be easily verified by means of (15), (16), 
and (18). These peculiar features may be removed by the following transforma- 
tion, which secures the same constant limits for both (20) and (21). 
For brevity put 



o 



_. , . ^, ^ cos / — cos /3 sin' ^/3 — sin^ \p 

Then, observme that = . , 1 ■21. . 

^ cos / — cos a sin'' \a — sm'' \p 

put in /i sin^/ = sin^ttsin?-!, 

and in I^ sin \p = sin \^ sin y^. 

_, . ,^ 2 si n ^a cos ^ dji 

These give dp = —j-, ^ . , 1 ^~2 — \ > 

^ ^ \/(i— sin' \a sin' ^i) 

, 2 sin -^i? cos 7-2 ^;-2 

'^^-l/Ci—sinH/Jsin'r^)' 

and the limits for both y^ and y^ are o and \Tt. Therefore 



i'^2sinii94( I I — ^!",?xsin'?-i 1 dy■^ 

/, =. r — 7^ ..r.^^ , ^ . (24) 

' J |/(i — sin'^asin^T-j) 



«</3 



. r 2 sin' ^^ cos' ra <^r2 (25) 

■ i«^[i-j5^sin'r2].l/(i-sin'i/3sin'r2)" 



o sin fa. 



«>/3 
[to be continued.] 



